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INTRODUCTION 


work  by  J.  Morrow  and  E.  S.  Rowland^  ^  aptly  described  the  sffoct 
of  residual  stress  in  fatigue  as  analogous  to  the  effect  of  mean  stress  and 
showed  that  the  effect  of  tensile  mean  stress  is  to  reduce  fatigue  life  while 
that  of  compressive  mean  stress  is  to  increase  fatigue  life.  It  has 
subsequently  been  incorrectly  implied  by  others  that  the  effect  of  residual 
stress  is  the  same  as  that  of  an  applied  mean  stress.  The  following 
relationships  show  that  they  are  analogous  but  not  the  same.  Considering  the 
mean  and  alternating  components  of  an  applied  cyclic  stress  gives: 

°max  =  °mean  +  °alt 


°min  °raean  °alt 


(1) 


However,  combining  a  residual  stress  with  an  applied  cyclic  stress  differs 


because  the  residual  stress  is  not  a  function  of  the  active  loading.  Thus 

°max  =  °residual  +  (°mean  +  °alt) 


°min  “  °residual  +  (°mean  "  °alt) 

l  ence  for  an  active  loading  from  zero  to  maximum,  as  in  a  pressure  vessel 
where  amean  is  equal  to  the  active  aalt : 


(2) 


°max  ~  °residual  +  2  oait 
°min  =  °residual 


(3) 


Morrow,  JoDean  and  Sinclair,  G.  M. ,  "An  Analysis  of  Cyclic  Stress  Behavior 
for  Conditions  of  Controlled  Strain,"  Theoretical  and  Applied  Mechanics 
Report  No.  543,  University  of  Illinois  Eng.  Bulletin,  August  1957. 

^Morrow,  JoDean  and  Millan,  J.  F.,  "Influence  of  Residual  Stress  on  the 
Fatigue  of  Steel,"  Society  of  Auto.  Eng.  TR-198  (July  1961),  SAE  J783,  Div.  4 
of  the  SAE  Iron  and  Steel  Technical  Committee. 

Rowland,  E.  S.,  "Effect  of  Residual  Stress  on  Fatigue,"  Proceedings  of  Tenth 
Sagamore  Conference,  Syracuse  University  Press,  1964,  pp.  229-244. 
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The  residual  stress  is  continually  present,  whether  or  not  the  service  load  is 
acting.  It  must  be  kept  in  mind  that  a  residual  stress  state  may  be  biaxial 
at  the  surface  of  a  body  and  may  be  triaxial  in  the  interior.  The  principal 
axes  of  this  state  may  or  may  not  coincide  with  the  axes  of  the  principal 
applied  stresses.  These  stationary  aspects  of  residual  stresses  make  it 
necessary  to  exercise  care  when  including  them  in  analytical  expressions  for 
fatigue  crack  growth.  This  report  is  intended  to  illustrate  a  way  that  the 
residual  stress  effect  may  be  examined  separately  from  mean-stress  or  R-ratio 
effects  in  fracture  mechanics  analysis. 

FRACTURE  MECHANICS  ANALYSIS 

Using  the  Paris-Erdogan^  fatigue  crack  propagation  rate  equation 

1 

da/dN  =  -  ( AK)m  (4) 

M 

one  can  include  the  effect  of  residual  stress  in  the  fatigue  life  estimate  for 

a  component  as  illustrated  in  the  following  analysis. 

The  crack  tip  stress  intensity  factor  K  resulting  from  an  applied  loading 

which  varies  from  zero  to  the  stress  level  a  is  given  by 

K  =  Ya/xa  (5) 

In  a  tension  member  a  =  P/A;  in  a  plate  or  beam  having  a  section  modulus  Z,  a 

is  the  applied  extreme  fiber  stress  Pe/Z,  where  Pe  is  the  bending  moment;  and 

’  2w2 

in  an  internally  pressurized  cylinder  of  diameter  ratio  w  =  OD/ID,  a  =  (-r — )p 

wz-l 


^Paris,  P.  C.,  "The  Fracture  Mechanics  Approach  to  Fatigue,"  Fatigue,  An 
Interdisciplinary  Approach,  Burke,  J.  J. ,  Reed,  N.  L. ,  Weiss,  V.  ,  Editors, 
Proceedings  Tenth  Sagamore  Army  Materials  Research  Conference,  Syracuse 
University  Press  (1964),  pp.  107-132. 
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from  the  hoop  stress  and  the  pressure  in  the  crack,  as  shown  by  Underwood.-* 
When  the  values  for  o  and  y  are  taken  to  be  independent  of  crack  depth  "a" 
these  expressions  apply  only  for  shallow  crack  depths. 


CRACK  SHAPE 

For  semi-elliptical  surface  cracks  we  define  a  shape  factor  as 

SF  =  --—  (6) 

1.12 

where  for  a  shallow  straight  fronted  crack  in  a  plate  under  tension  the  value 
of  y  is  1.12,  and  the  value  of  SF  is  unity.  The  values  for  y  for  curved 
fronted  semi-elliptical  surface  cracks  may  be  calculated  from  expressions 
developed  by  Newman  and  Raju^  for  surface  cracks  in  a  plate  under  tension  or 
bending  loads.  The  values  for  SF  for  a  given  aspect  ratio  and  crack  depth 
are  given  in  Table  I. 

As  the  crack  depth  increases  the  shape  factor  for  any  given  crack  shape 
increases  under  tension  loading  while  under  bending  the  shape  factor 
decreases.  Since  most  of  the  fatigue  life  of  a  component  is  expended  at  shal 
low  crack  depth  we  may  use  the  values  of  SF  for  a  «  B  as  constants  indepen¬ 
dent  of  crack  depth  for  a  good  first  approximation  in  either  tension  or  bend¬ 
ing  in  the  following  analysis.  Considering  SF  as  independent  of  the  crack 
depth  permits  a  relatively  simple  integration  of  da/dN.  More  precise  life 
estimates  may  be  made  by  taking  into  account  the  variation  of  SF  as  crack 
depth  increases. 


■^Underwood,  J.  H.  ,  ’’Stress  Intensity  Factors  For  Internally  Pressurized  Thick 
Wall  Cylinders,”  Stress  Analysis  and  Growth  of  Cracks,  Proceedings  of  1971 
National  Symposium  on  Fracture  Mechanics,  Part  I,  ASTM  STP  513,  American 
Society  for  Testing  and  Materials,  (1972),  pp.  59-70. 

^Newman,  J.  C.  and  Raju,  I.  S. ,  ’’Analysis  of  Surface  Cracks  in  Finite  Plates 
Under  Tension  or  Bending  Loads,”  NASA-TP-1578 ,  (1979). 
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TABLE  I.  SHAPE  FACTOR  ( SF)  FROM  NEWMAN  AND  RAJU  EQUATIONS6 


Crack  Front  Shape 

Aspect  Ratio 
a/ 2c 

Tension  or 
Bending  SF 
@  a  «  B 

@  a  = 

Tension  SF 

0.2W 

Bending  SF 

Straight 

0 

1.000 

1.000 

1.000 

Very  long  curve 

.005 

0.999 

1.144 

0.888 

Long  curve 

.10 

0.937 

0.993 

0.757 

Semi-elliptical 

.33 

0.719 

0.732 

0.543 

Semicircular 

.50 

0.586 

0.591 

0.432 

NOTE:  For  very  small  crack  depth,  (a  «  B) ,  the  value  of  SF  is  given  by 


1  1.13  -  0.09(a/c) 

SF  - - [ - -  ] 

1>13  /I~+  1.464(a/c)T-65 

for  either  tension  or  bending. 


INCLUDING  RESIDUAL  STRESS 

For  residual  stress  which  decreases  linearly  with  depth  from  the  value  0o 
at  the  surface  of  the  part  to  a  zero  value  at  the  depth  a 0  as  shown  in  Figure 
1 ,  the  crack  tip  stress  intensity  factor  is  given  by  Underwood  and  Throop^  as 

a  _ _ 

K  =  (1.12  -  0.68  — )  a0/n7  (7) 


^Newman,  J.  C.  and  Raju,  I.  S.,  "Analysis  of  Surface  Cracks  in  Finite  Plates 
Under  Tension  or  Bending  Loads,"  NASA-TP-1578 ,  (1979). 

^Underwood,  J.  H.  and  Throop,  J.  F.,  "Surface  Crack  K-Estimates  and  Fatigue 
Life  Calculations  in  Cannon  Tubes,"  Part-Through  Crack  Fatigue  Life 
Predictions,  ASTM  STP  687,  J.  B.  Chang,  Ed.,  American  Society  for  Testing  and 
Materials,  (1979),  pp.  195-210. 
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for  a  straight  fronted  crack.  Applying  a  shape  factor  for  a  serai-elliptical 
surface  crack  subjected  to  the  superposition  of  applied  stress  A a  and  residual 
stress  o0  gives 


°o  0.68 o0  a 

«max  -  SF  { 1.12(1  +  --)  -  ---- - }Aa/^  (8) 

Ao  (Aa)  aQ 


which  is  of  the  form 

Kmax  =  SF  {A  +  BajAa/rra  (9) 

If  there  is  no  residual  stress,  oQ  =  0,  A  =  1.12  and  B  =  0.  When  the 
residual  stress  is  compression,  0o  is  negative,  hence  A  is  smaller  than  1.12 

and  B  is  positive,  with  the  result  that  y^ax  is  smaller  than  when  there  is  no 

residual  stress.  Conversely,  when  the  residual  stress  is  tension  aQ  is 

positive,  hence  A  is  larger  than  1.12  and  B  is  negative,  with  the  result  that 
Kmax  is  larger  than  when  there  is  no  residual  stress. 


CRITICAL  CRACK  DEPTH 

The  critical  crack  depth,  ac,  for  a  given  applied  stress  range  Ao  can  be 
solved  from  the  expression  for 

1  Kic 

ac  {A  +  Bacp  *  -  [  7“~“  P  (10) 

7T  (SF)Ao 


Compression  residual  stress  increases  ac,  and  tension  residual  stress 
decreases  ac  relative  to  that  for  0o  =  0.  The  residual  stress  effect  is 
expressed  in  the  squared  term  (A  +  Bac) .  The  solution  of  equation  (10)  for  ac 
gives  three  roots  of  a  cubic  equation,  of  which  the  smaller  is  the  limiting 
crack  depth.  In  this  expression  A a  is  the  excursion  of  stress  from  zero  to 
its  maximum  value  in  the  load  cycle. 
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STRESS  INTENSITY  FACTOR  RANGE 


In  calculating  the  crack  propagation  rate  when  0o  is  compression  the 
stress  intensity  factor  range  is 


0q  0.68oo  a 

AK  =  (Kmax-0)  =  SF  { 1 . 1 2(  1  +  — ) - jAa/ira  (11) 

Ao  (Ao)  aQ 


or 


AK  =  SF{A  +  Ba  } Aa/ira  (12) 

Here  AK  is  the  same  as  Kmax  in  equations  (8)  and  (9)*  Since  o0  is  negative 
the  compressive  residual  stress  reduces  AK,  and  reduces  the  fatigue  crack 
propagation  rate  and  thereby  extends  the  fatigue  life  beyond  that  for  zero 
residual  stress.  If  Ao  is  not  great  enough  to  overcome  the  compressive 
residual  stress,  the  resulting  AK  will  be  negative,  the  crack  will  remain 
closed  at  maximum  load,  and  no  fatigue  crack  growth  should  take  place. 

When  the  residual  stress  oQ  is  tension,  i.e.,  positive,  the  stress 
intensity  factor  range  is 


AK  =  Kmax  -  Kfljin  (13) 

where  Kmax  is  given  by  equation  (8)  and  K^n  is  given  by  equation  (7) 
multiplied  by  the  shape  factor  SF.  The  result  is  that  for  tensile  residual 
stress 

AK  =  1.12(SF)Ao/iTa  (14) 


which  is  the  same  as  for  zero  residual  stress. 

Thus,  under  tensile  residual  stress  AK  depends  only  on  the  range  of 
applied  stress  and  the  crack  depth,  and  is  independent  of  the  tensile  residual 
stress.  Therefore  the  fatigue  crack  propagation  rate  is  relatively  unaffected 
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by  tensile  residual  stress.  Experiments  by  Underwood^  have  shown  that  tensile 
residual  stresses  produced  by  applying  compressive  overloads  do  not  greatly 
affect  the  measurable  crack  growth;  their  affects  are  apparently  on  the 
initiation  and  early  growth  of  the  crack.  However,  since  the  critical  crack 
depth  is  reduced  by  the  tensile  residual  stress  the  fatigue  life  is  reduced, 
but  the  effect  is  not  as  great  as  the  increase  in  life  that  would  be  caused  by 
an  equal  value  of  compressive  residual  stress.  Moreover,  tensile  residual 
stress  adds  to  the  applied  mean  stress  and  there  is  a  recognized  effect  of 
mean  stress  on  fatigue  crack  propagation  rate  which  is  not  treated  here  in 
this  analysis. 


S-N  EQUATION 

Combining  equation  (4)  and  equation  (12)  and  assuming  m  =  3  gives 


da/dN  =  -  (SF{A  +  BaUo/nl)3 
M 


(15) 


or 


Nf  M  1  M 

f  dN  - - -  J 

Ni  (SF/tT)3  ( A°) 3  a± 


da 

{a  4*  Ba  }3a3/2 


Integrating  as  suggested  by  Berge  and  Myhre 
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this  is  written  as 


N  =  (Nf  -  Ni) 


C  I 

(To)1 


(16) 


(17) 


^Underwood ,  J.  H.  and  Kapp,  J.  A.,  "Benefits  of  Overload  for  Fatigue  Cracking 
at  a  Notch,”  Fracture  frfechanics  Thirteenth  Conference,  ASTM  STP  743,  Richard 
Roberts,  Ed . ,  American  Society  for  Testing  and  Materials ,  1981,  pp.  48-62. 
^Berge,  S.  and  Myhre,  H. ,  "Fatigue  Strength  of  Misaligned  Cruciform  and  Butt 
Joints,”  IIW  Doc.  XIII  683-77  and  Norwegian  Maritime  Research  No.  1-1977, 
Norwegian  Inst,  of  Technology,  University  of  Trondheim,  Norway. 
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where 


®f  da 

I  -  /  - -  (18) 

ai  {A  +  Ba}3a3/2 

and 


M 

C  - - 

(  SF/tt)  3 

which  may  be  expressed  in  log-log  form  as 


(19) 


1  1 

Log  Aa  =  -  Log  N  +  -  (Log  C  +  Log  I) 


(20) 


Here  the  factor  -  -  is  the  initial  slope  of  the  log-log  plot  and  the 
1  3 


expression  +  —  (Log  C  +  Log  I)  is  the  intercept  at  N  =  1  on  the  graph  of  a 
versus  N.  (Note:  if  SF  is  not  independent  of  a,  the  terra  SF3  must  be 
included  inside  the  integral.  This  may  be  necessary  for  a  correct  solution 
for  deep  cracks.) 

When  there  is  no  residual  stress,  oQ  =  0,  and  A  =  1.12  and  B  =  0  in 

1 

equation  (18)  for  I.  Therefore  the  S-N  graph  has  a  uniform  slope  of  -  -  and 
an  intercept  of 


2  r  1  1  i 

-  (Log - +  Log  - - -  [  - - ]) 

3  (SF/?)3  (1-12) 3  —  — 


(20) 


1  M 

—  +  Log - -  [  - 

-x  3  n  nU  L 

a±  y  af 

When  the  residual  stress  is  compression  the  value  of  A  is  smaller  and  the 

value  of  B  is  larger,  and  the  limiting  crack  depth  ac  is  larger.  The 

resulting  value  of  1  is  increased  and  the  intercept  is  raised.  The  initial 

1 


slope  remains  the  same,  -  -,  but  for  decreasing  Aa  the  slope  may  change 
because  the  increase  of  critical  crack  depth  ac  may  increase  the  value  of  af. 
Both  effects  result  in  longer  life  at  any  Aa  when  compressive  residual  stress 
aQ  is  present. 
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When  the  residual  stress  is  tension  the  integration  of  equation  (4) 

1 

combined  with  equation  (14)  results  in  the  same  slope,  -  -  ,  as  for  zero 
residual  stress  but,  because  the  failure  crack  depth,  ac,  at  any  Aa  is  reduced 
by  the  tensile  residual  stress,  the  value  of  1  is  smaller.  Therefore  the 
intercept  of  the  S-N  graph  is  lowered.  The  net  effect  of  tensile  residual 
stress  is  a  shorter  fatigue  life  for  any  Ao  than  when  no  residual  stress  is 
present,  but  the  effect  is  predicted  to  be  less  than  for  equal  compressive 
residual  stress. 

These  effects  are  illustrated  on  the  following  schematic  log-log  S-N 
graph  and  in  the  following  numerical  example.  This  example  and  the  graph  were 
calculated,  using  the  equations  given  here,  expressly  to  compare  the  relative 
effects  of  compression  and  tension  residual  stresses  on  the  S-N  curve. 
Therefore  the  life  estimates  do  not  include  many  other  aspects  that  should  be 
considered  in  a  complete  fatigue  life  calculation. 

NUMERICAL  EXAMPLE 

Calculate  N  =  (Nf  -  N^)  for  Ao  =  100  Ksi  and  for  Aa  -  65  Ksi  in  bending, 

given  a^  =  0.002  inch;  a/2c  -  0.10;  SF  =  0.937;  and  Kqc  “  80  Ksi/in  ;  with  m  = 

3  and  -  =  3.4  x  10~10  for  martensitic  AISI  4340  steel  with  stress  in  Ksi  and 
M 

crack  depth  in  inches.  Compare  results  for  zero  residual  stress,  aQ  =  0,  and 
for  residual  stress  with  aQ  =  0.45  inch  and  tensile  aQ  =  +  40  Ksi  and 
compression  aQ  =  -40  Ksi.  The  results  are  given  in  Table  II,  and  the 
corresponding  S-N  curves  are  shown  on  the  graph  in  Figure  2. 
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TABLE  II.  CALCULATED  CRITICAL  CRACK  DEPTHS  AND  FATIGUE  LIVES 


Stress  Range 
Ksi 

Zero  Residual 
00  =  0 

Tensile  Residual 
o0  =  +  40 

Compressive  Residual 

0O  =  “  40 

Aa  =  100 

ac  =  0.185  in. 

ac  =  0.105  in. 

ac  =  0.318  in. 

Nf  =  4.34xl04 

Nf  =  4.19xl04 

Nf  =  1 . 56x10s 

Aa  =  65 

ac  *  0.438  in. 

ac  =  0.211  in. 

ac  =  1.177  in .  >  aQ 

*af  =  0.45  in.  =  aQ 

Nf  =  1 . 65xl05 

Nf  =  1 . 60xl05 

Nf  =  2.03xl06 

1  inch  =  25.400  mm;  1  Ksi  =  6.8948  MPa. 

*NOTE:  Since  Equations  (7),  (8),  and  (10)  only  apply  to  linear  variation 

of  residual  stress,  when  ac  calculated  from  Equation  (10)  is 
greater  than  aQ  for  a  compressive  residual  stress  the  value  of  af 
in  Equations  (11),  (16),  and  (18)  is  set  equal  to  aQ.  It  is 
thereby  assumed  that  the  number  of  cycles  to  grow  from  aQ  to  af  is 
negligible . 


S-N  GRAPH 


The  log-log  graph  in  Figure  2  shows  that  for  the  stress  range  between  65 
to  100  Ksi  (448  to  690  MPa)  the  tensile  residual  stress  causes  very  little 
reduction  of  fatigue  life,  but  the  compressive  residual  stress  causes  consid¬ 
erable  increase  in  fatigue  life  compared  to  the  zero  residual  stress  condi¬ 
tion.  Also  the  graph  for  compression  residual  stress  is  tending  to  curve 
toward  a  fatigue  limit  at  about  55  Ksi  (379  MPa)  stress  range.  In  real 
materials  these  tendencies  are  known  to  affect  both  the  crack  propagation 
threshold  and  the  fatigue  crack  rate  and  fatigue  life. 
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This  analysis  is  approximate  because  it  does  not  include  effects  of  the 
applied  mean  stress  or  the  changes  in  either  shape  factor  or  residual  stress 
distribution  as  functions  of  the  depth  of  crack*  However,  the  latter  two 
changes  are  nearly  negligible  in  the  very  early  portion  of  crack  growth, 
wherein  most  of  the  fatigue  life  is  spent,  and  the  solution  is  therefore 
meaningful • 

CONCLUSION 

Much  work  remains  to  be  done  in  developing  quantitative  relationships  for 
expressing  the  effects  of  residual  stresses  in  fatigue  for  purposes  of 
obtaining  accurate  fatigue  life  estimation  and  for  use  in  design  procedures. 
Many  aspects  of  residual  stress  effects  in  fatigue  behavior  have  been  merely 
mentioned  in  this  report  and  many  other  aspects  worthy  of  consideration  have 
been  omitted  entirely* 

The  relaxation  of  residual  stresses  during  the  service  life,  vibratory 
stress  relief,  effects  of  peak  overloads  in  service,  effects  of  prior  damage, 
residual  stress  effects  on  initiation  of  fatigue  damage  or  cracks,  contact 
fatigue,  and  heat  treatment  effects  are  subjects  that  have  been  suggested  for 
investigation*  Furthermore,  the  specification  of  methods  for  evaluating  the 
effects  of  residual  stress  on  fatigue  performance  and  for  controlling  the 
magnitude  of  residual  stresses  within  permissible  limits  need  to  be  addressed. 
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Residual  Stress 


Figure  1. 


The  idealized  variation  of  residual  stress  with 
depth  x  from  the  surface  of  a  component. 
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Figure  2.  A  fatigue  S-N  curve  on  a  log-log  plot  of  Act  versus  N  illustrating 
the  calculated  idealized-effects  of  a  tensile  residual  stress  and 
an  equal  compressive  residual  stress  compared  with  zero  residual  stress 
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